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Abstract. 

We study dynamics of entropy perturbations in the two-field assisted dark energy 
model. Based on the scenario of assisted dark energy, in which one scalar field is 
subdominant compared with the other in the early epoch, we show that the entropy 
perturbations in this two-field system tend to be constant on large scales in the early 
epoch and hence survive until the present era for a generic evolution of both fields 
during the radiation and matter eras. This behaviour of the entropy perturbations 
is preserved even when the fields are coupled via kinetic interaction. Since, for 
assisted dark energy, the subdominant field in the early epoch becomes dominant 
at late time, the entropy perturbations can significantly influence the dynamics of 
density perturbations in the universe. Assuming correlations between the entropy 
and curvature perturbations, the entropy perturbations can enhance the integrated 
Sachs- Wolfe (ISW) effect if the signs of the contributions from entropy perturbations 
and curvature perturbations are opposite after the matter era, otherwise the ISW 
contribution is suppressed. For canonical scalar field the effect of entropy perturbations 
on ISW effect is small because the initial value of the entropy perturbations estimated 
during inflation cannot be sufficiently large. However, in the case of k-essence, the 
initial value of the entropy perturbations can be large enough to affect the ISW effect 
to leave a significant imprint on the CMB power spectrum. 
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1. Introduction 

Observations, such as SNe la [I], WMAP |2j, SDSS [3], etc., indicate that the expansion 
of the universe is accelerating at late time jU [5]. Within the framework of general 
relativity, this acceleration of the universe may be explained by supposing that the 
energy density of the universe is dominated by an unknown form of energy, whose 
pressure is sufficiently negative, known as dark energy. It has been suggested that 
dark energy may be in the form of cosmo logical constant, exotic fluid, evolving scalar 
fields, etc. In general, scalar field models of dark energy, such as quintessence [6]- [8] 
or k-essence [9j [10], may be motivated by high energy physics theory. For example, 
quintessence may be a scalar field with the exponential potential V(4>) = Voe~ x ^, that 
arises in superstring and Kaluza-Klein theories [11]. Here, A is a constant parameter 
and the reduced Planck mass is set to be unity. 

For a single quintessence field with an exponential potential p^2l [13] . the evolution 
of the field in the early epoch will converge to the scaling solution for a wide range of 
initial conditions if A > a/3(1 + Wf), where wj is the equation of state parameter of the 
dominant fluid. When the scaling solution is reached, the energy density of the field 
will evolve in the same way as that of the dominant fluid [121 [13]. Hence, although the 
scaling behaviour may provide a solution to the cosmic coincidence problem, it cannot 
lead to an accelerated expansion of the universe. The accelerated expansion of the 
universe can be achieved if the field evolves towards the field dominated solution with 
A < y/2. Unfortunately, the range A > a/3(1 + Wf) does not overlap with A < y/2, so 
that the accelerating solution cannot be reached from the scaling solution. 

In general, there is no reason why the present accelerated expansion of the universe 
should be driven by a single scalar field. It has been shown that, for multiple scalar 
fields with exponential potentials, although the potential of each field is too steep to 
drive an accelerated expansion, the cumulative dynamics of these fields can lead to an 
accelerated expansion of the universe [H]. Based on this idea, multiple scalar fields with 
exponential potentials, with the slope of each potential in agreement with the prediction 
from particle physics, can drive an accelerated expansion of the universe. This idea has 
been originally proposed in the assisted inflationary scenario, and has also been used in 
dark energy model. This model of dark energy is assisted dark energy [151 CE] , and in 
this model, the scaling regime can be followed by accelerated epoch. 

Based on the motivation from high energy physics, it is possible that there exists 
an internal structure in the field space of multiple scalar fields, e.g., in DBI models 
[17] . Some features of this internal structure, which may be represented by Kinetic 
interactions between fields, have been studied using a simple phenomenological form of 
kinetic interaction in several dark energy models [18] - [21]. It has been shown that 
allowing kinetic interaction in assisted dark energy can have only minor effects on its 
background dynamics [T9] . 

In addition to the background dynamics, the features of dark energy also depend on 
the dynamics of its perturbations. In some cases, the perturbations in dark energy can 
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leave interesting imprints on CMB anisotropies. For example, in the case of single scalar 
field dark energy, the entropy (isocurvature) perturbations in dark energy can modify 
the ISW effect such that the CMB power spectrum is suppressed at low multipoles, if 
the entropy perturbations in dark energy and the curvature perturbations are initially 
anti-correlated and dark energy field is nearly frozen through out the evolution of the 
universe [22J - [24]. In this case, the dark energy field is required to be nearly frozen 
otherwise the entropy perturbation in dark energy damps away before it can affect the 
ISW effect. However, we will show that, for two-field dark energy models, it is possible 
for the entropy perturbations to survive until the present epoch for a generic evolution 
of both fields in some case. 

The dynamics of perturbations in multiple scalar field system have been studied 
by many authors [25] - [29]. In many cases, the dynamics of the field perturbations 
are conveniently described by decomposing the field perturbations into adiabatic and 
entropy modes [26] - [29] . 

In this work, we are interested in the behaviour of the perturbations, in particular 
the entropy perturbations, in the two-field assisted dark energy model. We use the 
formulation in [26l [27] to decompose the perturbations in dark energy into adiabatic and 
entropy modes. We also allow kinetic interaction between the fields in our consideration. 

This paper is organized as follows. In section 2, we present the phenomenological 
Lagrangian of assisted dark energy including the kinetic interaction between fields. In 
section 3, we first present the evolution equations for perturbations in general two-field 
dark energy model, and then apply them to the assisted dark energy model with mixed 
kinetic terms. We end with the conclusion in section 4. 

2. Assisted dark energy with mixed kinetic terms 

The general action for multiple scalar field dark energy with noncanonical kinetic terms 
can be written as [27] 



where index 1 — 1,2, ...N denotes the I-th scalar field, the reduced Planck mass is set 
to unity and 



where Gu((p K ) is the metric in field space and is used to raise or lower the field index /. 
The summation convention on both field and space-time indices is implicit. The inverse 
metric G IJ can be computed from the relation GjkG kj = 5j. In the following, we will 
simplify our consideration by supposing that Gu is constant. From ([[]), one can derive 
the energy- momentum tensor for dark energy [27], 




(1) 



X 



-GuV^Vcp 



(2) 



T „u = Pg ^ + Px Gjj V /yv V , 



(3) 
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where subscript ,x denotes derivative with respect to X, and the evolution equation for 
the scalar fields, 

GwV^VV)+P/ = 0, (4) 

where the subscript ,j denotes derivative with respect to cp 1 . In the following, we will 
use subscript ,a to denote derivative with respect to any variable A, and use subscript 
,/ where I is the field index to denote derivative with respect to field (p 1 . 

In a spatially flat FLRW universe, the evolution equation for scalar fields takes the 

form 

¥ + (sh + y^j¥~ jr GlJp J = ' ( 5 ) 

where a dot denotes derivative with respect to time t, cp 1 is the homogeneous part of (p 1 
and the Hubble parameter H is given by 

2 



H 2 



1^2 p» ^ 



where the index s = r,m,d represents radiation, matter and dark energy respectively. 
From 03]), the energy density of dark energy is given by 

p d = 2XP )X -P. (7) 

Multiplying (jSJ) by (pi, we obtain 

I + 6 ^I + ^ = 0, (8) 

Pd,X 

where c 2 = Fx / p d ,x is the effective sound speed of dark energy. 

The important ingredient in assisted dark energy model is that the dark energy 
scales as the dominant fluid, i.e., undergoes the scaling regime, in the early epoch. This 
implies that the equation of state parameter Wd = P/pd of dark energy is constant, and 
therefore P oc pd in the scaling regime. From (jBJ), one sees that pd oc P oc H 2 oc t~ 2 . 
Hence, in the simplest case where Fx is constant during the scaling regime, it follows 
from (jTJ) that X oc t~ 2 , and consequently <p ! oc \jt and (p 1 oc mi, where (p 1 are the I-th 
fields that contribute to the dynamics of dark energy during the scaling regime. If we 
further suppose that the sound speed c e is also constant during the scaling regime, (jHJ) 
shows that Pd.i^p 1 oc t~ 3 , which is possible when the (//-dependent parts of P depend on 

Q=fc2cje*>A , (9) 

where cj and Aj are constant. During the scaling regime, Fx is supposed to be constant, 
so that we can write P oc Xg(X,Q), where g(X,Q) and Xgx(X,Q) are constant 
during scaling regime. A possible form of g(X,Q) is g(X,Q) = g(Y), where g(Y) is a 
polynomial function of Y, and Y is the combination of X and ([9]) such that it is constant 
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during scaling regime. From the above consideration, one sees that dark energy can scale 
as the dominant fluid in the early epoch if its Lagrangian takes the form 

P(X,<p I )=Xg(Y), (10) 

where, in general, Y is the sum of the following expressions 

X 



Y a = x[J2cje^ J ), Y b 



-\j<pJ 



where \j(p J > 0. 

Let us now compute Q d during the radiation and matter eras, and Wd during the 
dark energy era. Substituting ( ITU]) into (JSJ), we get 

- 6HP X X = Xp d , x + p dJ ^ = X (2X 2 g, xx + 5Xg, x + g) + X (2Xg, XI + 9j ) ft , (12) 

where the index / runs over the fields that contribute to the dynamics of dark energy. 
It can be seen that, during radiation and matter domination, the above equation can 
be satisfied if (j) 1 = (2/Aj)lnt+ constant and hence (p 1 = (2/A/t). Using this expression 
for cj) 1 , ( ITU]) and pip , one can show that 

Xg tX = -j9,i¥i and x 9,xx = l -^9,xi¥ ~ 9,x ■ (13) 
Substituting the above equations into (TT2]) . we get 

where 



^ = ^ = 3(l + ^/)^f, (14) 



l^ = G ir — -. (15) 
A eff A/A// 

In the above equation, the index I' also runs over the fields that contribute to the 
dynamics of dark energy. Similarly, one can show that Wd during dark energy domination 
is given by 

^ = -1 + 7^- (16) 

It follows from (IT5]) and (ITUj) that for suitable choices of G///, A e ff decreases and hence 
Wd gets closer to —1 when more fields contribute to the dynamics of dark energy. This 
shows that the multi-field dark energy described by (fTU]) has assisted behaviour. Note 
that the form of Lagrangian in ( flUi) does not coincide with that given in [16]. This 
is because the internal structure in field space as in (J2]) can give rise to the coupling 

terms between the I-th field and the other fields such as 7 </>/', [&J e A/< ^ (/' ^ I) 
etc.. These terms are absent in the Lagrangian in [16J. However, according to the above 
consideration, the scaling solution and assisted behaviour exist even when such terms 
appear in the Lagrangian. In general, the form of Lagrangian (fTU]) cannot be reduced 
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to that given in [TH] by choices in the forms of Gu, g{Y) or Y. Nevertheless, in some 
cases, the both Lagrangians can become identical (see e.g. (fl9l ). 

In the scenario for assisted dark energy, A e fr must be significantly larger than unity 
in the early epoch because during radiation era is required to be smaller than 0.045 
according to BBN constraints [30J. For accelerated expansion of the universe, A e g has 
to decreases at late time. This is possible when some of the subdominant fields, that 
had no contributions to the dynamics of dark energy in the early epoch, can start to 
contribute to the dynamics of dark energy at late time. 

We now show that although (fTUl) has assisted behaviour, not all the possible forms 
of Y presented above can realize the scenario for assisted dark energy. It can be shown 
that ([5]) can be written as 

if = - 3<^' - i||Uv* - + • < 17 > 

= - am 1 - - \tt/ p * + F ~ J TT- ■ (18) 

Here, (f) 1 are the I-th dominant fields that contribute to the dynamics of dark energy 
in the early epoch, and <p J are the J-th subdominant fields in the early epoch. For 
simplicity, we set G IJ = 5 IJ , Substituting (ITUi) into i fPTj) and ( |T8l) . one can see that, 
for all possible form of Y, all terms on the RHS of (ITT]) are proportional to t -2 
if (p 1 are in the scaling regime. Nevertheless, for some forms of g(Y) and Y, e.g., 
g(Y) = — 1 + Y = — 1 + Y a , the fourth term on the RHS of (fT8l) decreases with time 
much faster the the remaining terms and the fourth term on the RHS of ( jTTl) . We will 
consider this case below. Since J are subdominant, |0 J | <C |0 7 |. According to the 
scenario for assisted dark energy, the expansion of the universe will start to accelerate 
when some of J begin to be dominant. Hence, in order to drive a present accelerated 
expansion of the universe, \<fi J \ should be comparable with |0 7 | about the present epoch. 
It follows from flI7j and flE]) that the ratio \<j) J \ / Ifi 1 \ can increase if <p J /cj) J are less 
negative than (p 1 /^ 1 . For convenience, let the forth term on the RHS of (|T7|) and (TT8l) 
be denoted by Fj and Fj respectively. In the case where <fi J > initially, Fi and Fj 
are possitive while the remaining terms in their equations are negative. This implies 
that <p J / cj) J will be less negative than <f) ! /(p 1 if the initial conditions for <fr J are chosen 
such that Fj/(f) J > Fi/cf) 1 . Since Fj/<j) J decrease with time much faster than Fi/ip 1 and 
Fj decrease faster than the remaining terms in their equation, the ratio |0 J |/|0 7 | will 
increase only for a short period of time if the initial values of Fj / <p J are not significantly 
larger than Fj/^ 1 . If the initial values of Fj/<p J increase, the rate of growth of l^l/l^l 
will increase. Nevertheless, the period of increasing |0 J |/|0 7 | does not expand so much. 
As a result, J will start to be dominant in the early time when the initial value of F 



increase. Using similar consideration, one will see that <p J also cannot be dominant at 
the right time for the case where <p J < initially. 

In the following, we restrict ourselves to the case of two fields. In the early epoch, the 
evolution of entropy perturbations in this dark energy model can be studied analytically 
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using (TTOT) . However, numerical integration is required to study the effects of entropy 
perturbations on CMB power spectrum, so that we have to specify the form of P(X, <f) ! ). 
We consider both the cases of quintessence and k-essence whose forms of P(X, (J) 1 ) are 
presented in the following sections. 

2.1. The quintessence with exponential potential 

If we set g(Y) = 1 - Y = 1 - Y c , (HDJ) becomes 

P = X - e' Xl(t>1 - e~ x ^ 2 . (19) 
Here, we have set c\ — ci = 1 for convenience. If we choose the form of Gu as 

Gu = 8u + (5 n $j2 + SnSji) a , (20) 

where \a\ < 1 is constant, the Lagrangian ( [19]) becomes that considered in [19] and 
[T6] when a = 0. We will use this form of Gu in our numerical integration for both 
quintessence and k-essence models. 

2.2. The k-essence model 

We now set g(Y) = -1 + Y = -1 + Y a - M 4 Y C , where M is constant, so that (TTOl) 
becomes 

P = -X + X 2 (W 1 + c 2 e A2 * 2 ) - M 4 (V Al< ^ + c 2 e~ x ^ , (21) 

where we have set c\ in Y a and Y c equal to unity. The above Lagrangian has the form 
of the Lagrangian for dilatonic ghost condensate model studied in [31J. According to 
the previous consideration, g(Y) = — 1 + Y a cannot realize the scenario for assisted dark 
energy because the fourth term on the RHS of (Tl8|) decreases too rapidly compared 
with the others. Nevertheless, due to the term M 4 Y C in the expression for g(Y) in this 
section, the fourth term on the RHS of ffT8|) does not decrease much faster than the 
remaining terms. As a result, 2 can contribute to the dynamics of dark energy at late 
time, i.e., the assisted scenario can be realized. 

The value of C2/M 4 that can realize the scenario for assisted dark energy, depends 
on the value of M 4 , the starting time for the numerical integration and the contribution 
from (f 2 to the dynamics of dark energy at initial stage, i.e., the initial value of 2 /0 4 
and e A2( ^ 2_Al9il . In our consideration, we set M 4 = 1 and C2 = 10~ 40 . 

3. The perturbations in assisted dark energy 

We will mainly perform our calculations in the longitudinal gauge [32] , in which the line 
element is 

ds 2 = -(1 + 2^)dt 2 + a(t) 2 (l - 2<S>)5 ij dx i dx j = g^dx^dx" . (22) 

The evolution equation for the field perturbations dcj) 1 = y? 7 (£, x) — 7 (t) can be obtained 
by expanding (TJJ around the FLRW background. From (TJJ we write 

GWV M (PxV„0 J ) + Pj = . (23) 
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Expanding the field variable around its homogeneous background and substituting g^ v 
from (122]) into the above equation, we obtain the following linearized equation for Sep 1 
in Fourier space: 



P 



x 



P.x I a 



+ -I U^SP, X 
a 6 V 



= P^¥ (3$ + ^) + 2G IJ PjV + G IJ 5P )J , (24) 

where SP^x and SPj are the perturbations in P x and Pj respectively. In the single 
field case this equation is in agreement with [33j [31]. The evolution of the metric 
perturbations is governed by the perturbed Einstein equation [33] . 

3#(6 + #*) + J$ = -\J2 6 P°i ( 25 ) 

s 

$ + = X -ak- 1 J2 C 1 + ^) P-«»> > ( 26 ) 
= a 2 A;" 2 ]T ( Ps n s ) , (27) 

s 

where Sp s , v s and IT are the density perturbations, velocity perturbations and 
anisotropic perturbations in specie s respectively. For dark energy, ILj = and 

Sp d = 25 (XP X ) -SP = P X 5X + 2XSP, X - Pj6(f> J = p d , x SX + p dyJ 8<\> 3 , (28) 

(1 + w d )p d v d = p d u d = P^koT 1 ^ '84>j , (29) 

where 

8X =0 J 50 J -2X^, (30) 

SP = P X 5X + P yJ 5(j) J , (31) 

SPx = PxxSX + P X jS<p J ■ (32) 

Using f |24l) . f |28l) and ( l29l) . we obtain the evolution equation for S d = Sp d /p d and u d , 

/ftp \ . u 

5 d = -3H — -w d )5 d + 3(l+w d )$--u d , (33) 

kSP k 

u d = H (3w d - 1) u d + + {l + w d )-V. (34) 

a p d a 

As expected, these equations are the evolution equations for the perturbations in a 
perfect fluid whose pressure perturbations are SP. These equations can be closed if the 
relation between SP and Sp d is specified. For the the Lagrangian of the form P(X, if 1 ), 
it is convenient to compute this relation by decomposing the field perturbations into 
adiabatic and entropy modes. 

3.1. Decomposition into adiabatic and entropy modes 



The field perturbations can be decomposed into adiabatic and entropy modes using a 
set of suitable orthonormal bases in the field space. According to [261 EZ], we choose 



Dynamics of entropy perturbations in assisted dark energy with mixed kinetic terms 9 
the first basis as 



(35) 



V2X 

The other basis can be specified by solving the equations 

s 7 G 7 js J = s/s 7 = 1, and e'G/js 7 = e 7 s/ = . (36) 
Using these bases, the field perturbations can be decomposed as 

8(f) 1 = Sae 1 + SSs 1 . (37) 

From the orthogonality of these bases, the adiabatic mode 5a and entropy mode SS can 
be written in terms of the field perturbations as 

5a = e I G IJ 5(j) J = e J 5(j) J JJ and 5S = s,j5(f) J . (38) 

If this decomposition is applied to the background field, the perturbations 5(f) 1 should 
refer to a variation in time. Hence, we can write 

a = e J( p J = V2X, and S = sj0 J = . (39) 

This shows that the background entropy field is constant, which is in agreement with 
the fact that the "relative" entropy perturbations 5S is gauge invariant. 
Multiplying (jSJ) by ej, we obtain the evolution equation for a 

a = -[3H+^-)a + ^ = -3c 2 e Ha - B&L , (40) 
\ p ,x ) P,x Pd,X 

where P CT = P/e 7 . On the other hand, if we multiply (JS]) by s/, we will get a useful 
relation 

s 1 = -J^e 1 = -Oe 1 , (41) 

where Pg = Pjs 1 . Using the relation ejs J = 0, we obtain 

e 7 = ^s 7 . (42) 

The above equations can be used to compute the relation between 5P and 5pd- The 
pressure perturbations 5P in fl3Tl) can be written in terms of 5a and 5S as 

SP = P X 5X + P t Ja + P tS 5S . (43) 

If we also write Spa = pd,x5X + pd, a 5a + Pd,s5S and use the relation P x = c\p<i,x-> we 
will get 

5P = c\5 P d + (P CT - c\p d ^ 5a + (P s - c 2 ePdjS ) 5S . (44) 

Since S = 0, P CT — c\pd,a = (P — c 2 e Pd^ /a. Using the definition of the adiabatic sound 

speed c\ = P/pd, and the conservation equation for p d , pd = — 3H (1 + wj) Pd, the 
second term on the RHS of the above equation becomes 

3H (c 2 e - c 2 a ) P x a5a = 3H^ [c\ - c 2 a ) p d u d . (45) 
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Here, we have written u d in terms of the adiabatic field as 

p d u d = P } xka~ l 4> J 8(f)j = P^ x ka~ 1 &5o- . (46) 

From (J7J), we can write the third term on the RHS of (jSj) in terms of the pressure only. 
Finally, we have 

SP = c 2 e 5 Pd + 3^H (c 2 e - c 2 a ) Pd u d + [P s (1 + cl) - 2c 2 e XP xs ] SS . (47) 

If SS = 0, the form of this relation will be the same as the one for single scalar field 
derived in [22]. This relation implies that ( 1331) and ( |34|) can be closed if the evolution 
of S S is known. 

In order to derive the evolution equation for SS, the following relations [27] are 
required, 

Ps = ~ ?^ + P,*so- + P,xsv° , (48) 
oP,x 

P„s = Pue I s J J and P xs = P XJ s 3 . (49) 
The evolution equation for SS is computed by multiplying (1241) with sj. The result is 

P >x \ xc , f k2 , ..2 , - 2 



where 



*S i I oH + ^J 5S+ ^- +l4 + -jSS= -fip-Sp*, , (50) 

U 2 = _^S5 1 P2 S 2 P,X S P,S (51) 

P tX 2ciXP 2 x +2 P 2 X ' (51j 

S = "4- [(1 + C e 2 )P,5 - C^, XS( T 2 ] , (52) 

and is the density perturbations of dark energy in the dark energy rest frame, given 
by 

S Pdv = 6p d +3H^ Pd u d = p dtX (&6a - - a 5a) - P s + lj - 2XP XS 6S .(53) 
In the second equality, the expression of Sp d in terms of the adiabatic and entropy fields, 
S Pd = Pd ,xbX + p d ,j5<p J = Pd x (&5a - 2X^ + Pd< Ja - P s + X ) " 2XP ,xs] SS . 



has been used. For the early epoch, it is possible to solve (133[). ( l34j) and (j50|) . 
which completely describe the dynamics of perturbations in dark energy, analytically 
Nevertheless, numerical integration is required to study the behavior of perturbations in 
dark energy about the present epoch. In order to solve these equations numerically, the 
initial conditions for S d , u d , 5S and SS in the early epoch must be specified. Since the 
contribution from entropy perturbations SS is small in the early era, it is possible to use 
the following adiabatic initial conditions for Sd and u d if w d is approximately constant 
during the initial stage. 

and — — — = v r , (54) 



1 + w d 1 + w r l + w d 
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where 5 r = 5p r /p r and v r is the density contrast and velocity perturbations of radiation 
respectively. Supposing that dark energy is a fundamental field which has already existed 
since inflationary era, the initial magnitude of 5S may be estimated by studying the 
dynamics of perturbations in dark energy during inflation. 



3.2. The initial conditions from inflation 

Here, we would like to roughly estimate the initial magnitude of SS by studying the 
evolution of perturbations in dark energy during inflation. In the following calculation, 
we assume for simplicity that inflation is driven by a single canonical inflaton field, 
the energy density of dark energy is sub dominant, dark energy and inflaton are 
weakly coupled via gravity only. Since it is convenient to study the dynamics of 
perturbations during inflation using spatially flat gauge, we write the evolution equation 
for perturbations (12^1) in spatially flat gauge using the relations [36J, 137] 



flat 



^ + {h) ' Q i = W** = s <l>L s + <r 1 jr- ( 55 ) 

where subscripts "flat" and "long" denote spatially flat and longitudinal gauge 
respectively. After performing a transformation , the metric perturbations \I/fl a t can 
be eliminated using the relation [26] 

1 



flat 



2H 



(56) 



where ip is the background inflaton field and Qf is the perturbed inflaton field in spatially 
flat gauge. It is not hard to show that (HM in spatially flat gauge is compatible with 
the action 



S (2) = i 

2 



dtd 3 ka 3 



PxGrj + P 



xx9i<Pj 



Q'Q" 



Pxk' 



+PuQ I Q J - C K Q f Q K - D K QfQ K + 2P XJ <j )I Q J Q I \ , (57) 

whose form is similar to that in [27], except for the third, fourth and fifth terms in the 
action. This is because the metric perturbations \I/fl at in (156]) are sourced by the inflaton 
field not by Q 1 in our case. Thus, to simplify our task, we follow the calculation in [27]. 
Here, inflaton and dark energy are coupled via the perturbed metric, and the coupling 
coefficients are given by 



C 



2X 



K 



H 

The action (15 



2 



P,XK^ 



H 



Px[l + 



A 



D 



P 



x 



K 



H 



M ■ (58) 



can be written in terms of the adiabatic and entropy fields as 



dtd 3 ka 3 



P 



.x 



Pxk 2 
a 2 



QmQ 



N 



(59) 



+ P,mnQ m Q n ~ C M QfQ - DmQ/Q + 2a (Q a - 9Q S ) P^uQ 1 
where indices M, N run over a, S. Using the new fields which are canonically normalized, 



Qa = PQc 



v s 



a V^x Qs = iQs 



(60) 
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together with conformal time r = J da /a, the second-order action (1591) can be rewritten 
in the form 



drd 3 k [v' 2 
+0„ 



,'2 



- 2Zv' a v s - k 2 c 2 e v 2 - k v s (61) 
Vt ss v 2 s + 2tt aS v a v s + C a Q f v a + C s QfV S + D a Q' f v a ] , 



where a prime denotes derivative with respect to conformal time, and 
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The evolution equation for entropy mode can be derived from ( I6T1) as 



7" 



7 



F m ° 1 < 



a 2 a 



2c P 



H 



iQj ■ (63) 



In order to solve this equation, the terms that are propportional to £ and //| are supposed 
to be negligible, i.e., |£| <C aH and /i| <C if 2 , which is possible when one dark energy 
field is subdominant compared with the other. We will discuss this issue in the next 
section. According to the scenario for assisted dark energy, the field, which undergoes 
scaling regime during radiation era, has larger contribution to the energy density as 
compared to the other before it enters the scaling regime. Since, before the scaling 
regime starts, the subdominant field tends to be nearly frozen [TO] while the dominant 
one can undergo a long period of kinetic regime (see figure (JTJ)), the contribution of the 
dominant field to the energy density of dark energy does not tend to decrease when 
looking back in time. Therefore, it is possible that the field that is dominant during 
radiation era is also dominant during inflation. If |£| <C aH, fi 2 s <C H 2 and the dark 
energy fields slowly evolve, i.e., Fx and c e evolve slower than the scale factor, the 
approximate solution of fl63|) is 



v s {k) 



-ikr 



2k 



i 

kr 



For a slow rolling inflaton, it can be shown that [37] 

ay 2 k 

on large scales, so that 



(64) 



(65) 
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where 1Z is the curvature perturbations on the comoving hypersurfaces, e is the slow roll 
parameter of inflaton, R is the relative amplitude of the tensor to scalar perturbations, 

and || ■ • • || denotes yj ((• • -) 2 ) and | • • • | when Q s and Qf are uncorrelated and fully 
correlated respectively. Here, (• ■ ■) refers to the ensemble average. In the second equality, 
we have used the fact that 5S is gauge invariant. It follows from this equation that the 
magnitude of SS can be large if \P t x\ is small, which is possible when dark energy slowly 
evolves in many k-essence models, e.g., in the case of Lagrangian in ( )2T1) . Px — > when 
Wd — y — 1 if M 4 = 1. This result can be used to estimate the initial magnitude of 
SS during the radiation era, because, as we will see in the next section, SS tends to 
be constant for generic evolution of dark energy as long as one of the scalar fields is 
sub dominant. 

Since we use the Lagrangian of the form P(X, (p 1 ) for dark energy in the 
above calculation, the adiabatic and entropy mode of the dark energy perturbations 
propagates with the speed of sound c e and speed of light respectively. However, in 
more general case in which the Lagrangian takes the form P(X IJ , ip K ) [29] where 
X IJ = — (l/2)V Ai y3 / V A V J , the propagation speed of adiabatic and entropy modes can 
be different from the ones considered here. For example, in the case of DBI Lagrangian 
[28| 129] . both adiabatic and entropy modes propagate with the speed of sound, so that 
the gradient term in the second order perturbed action is multiplied by sound speed. 
In this case, the amplitude of the entropy field SS can be amplified compared with the 
adiabatic field if |c e | < 1. However, in our case, it follows from (|66|) that the magnitude 
of SS can be enhanced compared with the curvature perturbation when \Px\ < 1, 
because, in order to obtain the canonically normalized field variable, we rescale the 
entropy field SS by the scale factor and square root of is kinetic coefficient yTy, and 
the curvature perturbation is generated by canonical inflaton field. 

In general, it is not convenient to use the value of Q a during inflation to define 
the initial conditions for 5a in the radiation era, because the evolution of Q a depends 
on the evolution of dark energy after inflation. Hence, we define the initial conditions 
for the 5a and also 5a through the adiabatic conditions for the perturbations in the 
total energy density and velocity given in (15^|) . However, strictly speaking, (15^|) can 
be justified if the curvature perturbation and the perturbations in dark energy have 
correlation. Nevertheless, in the above calculation, the perturbations in dark energy and 
curvature perturbation are uncorrelated because we suppose that |£| <C aH, /z| <C H 2 
and the inflaton slowly evolves. Several models for generating correlated dark energy 
perturbations have been proposed [221 123]. To simplify our task, instead of considering 
these models in detail, we use the model in [23] to estimate how (16 6 p is modified when 
the perturbations in dark energy and curvature perturbation are correlated. Based 
on the model in [23], the correlated dark energy perturbations can be generated if 
the curvaton field is coupled with the dark energy fields such as gi^C where ( is 
the curvaton field and gj is the coupling coefficient. In the case where |£| <C aH, 
/i| <C H 2 and Px together with c e evolves slower than the scale factor, the dynamics 
of Qs and Q a should not be much different from the dynamics of perturbations in dark 
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energy field in [23J. Hence, comparing the result in [23] with (jHSJ), one expects that 
(166]) are modified only by multiplication factor whose magnitude is smaller than unity, 
such as \8S\ < y/0.125R\TZ\/ \fP~x-, when the entropy perturbations in dark energy and 
curvature perturbation are correlated. Thus, in the following consideration, we assume 
that the perturbations in dark energy and curvature perturbation are fully correlated 
and use the RHS of the last equality in f )66|) to give a rough estimate of the upper bound 
of \SS\. 

3. 3. Evolution of the entropy perturbations during the radiation and matter epochs 

In order to study the evolution of the entropy perturbations, the form of s 7 should be 
known. In the case of two dimensional field space, applying the constraint equations for 
basis s 1 in (1361) to a symmetric Gjj yields 

b1= C 22 e 2 + G_ 12 e^ ^ g2 = ± Gne^+£^ 
VdetG VdetG 
where detG > 0. Obviously, there are two possible choices of s 7 in the two 
dimensional field space. These choices correspond to two s 1 pointing in opposite 
directions. Nevertheless, these different choices of s 1 do not lead to different dynamics 
of perturbations in dark energy. This is because SS and the derivative d n /dS n , where 
n is an odd integer, change sign for these different choices of s 1 , so that the evolution 
equation for SS in ( 1501) . the relation between SP and Spa in ( )47|) and the expression for 
Spd in terms of adiabatic and entropy modes are the same for these different choices. 
However, it follows from (1661) that these different choices of s 7 imply a change in the 
relative sign between SS and the curvature perturbations. Nevertheless, this relative 
sign cannot be computed explicitly from the choices of s 1 . In our consideration, we 
choose s 1 = — (G 22 e 2 + Gi 2 e 1 )/y/detG and s 2 = (Gne 1 + Gi 2 e 2 )/A/det G, and set the 
ratio SS/1Z to be negative initially, i.e., SS and the curvature perturbation are anti- 
correlated. We set SS/1Z < initially because it leads to the suppression of the CMB 
power spectrum at low multipoles for our choice of Lagrangians (|T9|) and (|2ip in our 
numerical integration. However, in general, initial SS/1Z < does not necessarily lead 
to the suppression of the CMB power spectrum at low multipoles. 

In the radiation and matter epochs, when the field (f 2 is subdominant compared 
with 1 , the magnitude of s 1 can be estimated by solving fl5j. Rewriting fl5j as 

^(a 3 P x /) =a 3 G /J Pj. (68) 



(69) 



dt 

When <j) 2 is subdominant, \P t i\ 3> |-P 2 |, and hence 

where C\ and c 2 are the constants of integration. This result implies that 



' + ZTf^ c 2 + 7^ c i • ( 7 °) 



G 22 a 3 P ' x V G22 
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Here, we have used 
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12 



G 



12 



detG 

Substituting (1701) into (j2J), we get 



and G 
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detG 
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12 



G 



22 



-i 2 



(71) 



(72) 



In the case of our interest where — 1 < Wa < 1, the first term on the RHS of this equation 
will not decrease faster than the second term, so that the second term will be negligible 
compared with the first term as the universe evolves. Thus, we ignore this term in the 
following calculation. Substituting ( 1701) into ( 1671) . we obtain 



< s z 



VdetG , 
-G—° 



22 



where 



G 



22 



detG 



(73) 



Here, the expression for | 1 is computed using (172"]) and (155]) . In the following 



calculations, we suppose that |s 2 



0(1). 



Let us now come to the case of assisted dark energy whose Lagrangian is given by 
( TlP|) . Using this Lagrangian, we estimate the magnitudes of //| and H which correspond 
to the amplitude of the "effective mass" of SS and the coupling between 5S and 5prf. 
From ( I13p . one can derive the following relations in the case where <p 2 is subdominant: 



\\x\Xg )X ^ 0,1 , and 
Using these relations, |Pi| ^> |P 2 | and 
can be seen that 



\9,u\ 



(74) 



3J), if Ai is not very much larger than unity, it 



X\P >xs \ = X\(g + Xg 



x. 
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Ai 



,ss\ 



IsVPu 

4ol 



2s 1 s 2 Pi 2 + s 2 s 2 P 22 | 



IP.11 5. 



11-DSS 



~ |AiPiPs5 

IP55I 



(75) 

~ |AiP j(T P S s|,(76) 



<C 1 when 



is 



where B ss = sV + 2sVPi 2 /Pn + s 2 s 2 P 22 /Pn so that 
subdominant. The relation (175]) is valid when Pxs 7^ 0. From Pi ^> Ppand (173]) . 
it is easy to see that 



|P,S| = IPi^+PoS 2 



IP.Pc 



(77) 



1^,1^1 

where B$ = s 1 + s 2 P 2 /Pi so that |Ps| <C 1 when 2 is subdominant. Using ( 174]) . (jTOl) 
and (1 + w d )p d = 2XP tX , we get 



Pd 

Hence, 



\iXg. 



x 



(1 - w d ) 



Ai 



and 



Pr 



Ai 



1 - Wd 



(7* 



2X 2 ^ X + X(7 2X 9)X + g ^ 2 ' XP X "^1 + ^ 

jPsl/pd ~ Ai(l — u>d)|Ps|/2 <C 1, implying that the contribution of to SP is 
negligible in the early epoch. In the early epoch, when the contribution from 2 to p d 



is significantly small such that \B>ss\ 
(T77|) and d78D that 
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,5 



crPx 
P,ss 



Pdp, 
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< 
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P.P.c 



<C 1 + w d and IP5I <C 1 + w d , it follows from ( 1761) 
A1P5 1 
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(79) 
(80) 
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ia 3 P tX SS) wO. 



In the above calculations, we have used \&\ < yfpd instead of \&\ ~ ^fp~d-, because \&\ can 
be much smaller than y/p2 when Wd — > — 1 in some cases, e.g., in the case of quintessence 
models where P = X-V. Using (jZHJ), flj2|), (EO]) and ^ < # together with c 2 ~ 0(1), 
it can be shown that |S| <C H and /i| <C if 2 . 

When |5| <C H and /i| -C if 2 , the evolution equation for SS ( 1501) on super horizon 
scales becomes 

d_ 
dt 

The solution of this equation is 

• di f didt SS , . 

Sb = - - , Sb = — - h const. ~ — — + const , (82) 

a 6 P yX J a A P,x H 

where d\ is the constant of integration. To evaluate the integral in the above expression 

for 5b, Fx is assumed to be approximately constant. Since, according to the above 

equations, \5S/H\ decreases with time, and \5S/H\ < \SS\ on large scales during 

inflation because the dark energy fields slowly evolve, we get SS ~ constant on large 

scales after the inflationary era. 

To check the validity of the above consideration, we solve ([5]), ( 1331) . (1341) and 

(150]) numerically using the modified version of CMBEASY [38]. In the following, 

the cosmological parameters for the background are chosen such that = 0.044, 

ftW M = Q.216, = 0.74, h = 0.72. In our numerical integration, the form of Gu is 

chosen according to (1201) . so that 

and 2 - + ^ (83) 



vT^ 2 vT^" 



a 



We first consider the case where a = using the quintessence model whose 
Lagrangian is given by f|T9|) . The evolutions of Wd = P/pd and w 2 , which is the equation 
of state parameter of subdominant field 2 , for various initial conditions are plotted in 
figure ([1]). There are four types of the initial conditions in figure ([I]). These four types 
correspond to the four possible combinations of the thick and thin lines in this figure, 
e.g., the thick solid together with thin long dashed lines correspond to the case where 
dark energy is initially in the potential regime Wd ~ —1, while 2 is in the kinetic regime 
W2 ~ 1 initially. All of these initial conditions, except the case where Wd ~ — 1 while 
u>2 ~ 1 initially, lead to the same evolution of 5S in the early epoch, as shown by thin 
solid line in figure ([3]), This evolution of 5S is also similar to that in the case of a ^ 
represented by thick solid line in the same figure. In the exceptional case where Wd ~ — 1 
while u>2 ~ 1 initially, the kinetic energy of 2 can be larger than that of 1 initially 
although 2 gives a negligible contribution to the energy density pd, so that |e 2 | > | e 1 1 
at the initial stage, and therefore the above analytical consideration is not applicable. 
The evolution of 5S for this case is represented by thin long dashed line in figure ([3]). 
This figure shows that 5S decreases until w 2 is close to —1, i.e., | e 1 1 becomes larger than 
|e 2 |. Due to the source term on the RHS of ( |50l) . 5S does not decrease to zero. Thus, 
this case does not imply the vanishing of SS. 
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We now consider the case where a^O using the k-essence model whose Lagrangian 
is given by (|2ip . The initial conditions for k-essence in figure (T5]) are chosen such that 
either <fi 2 = ia^ 1 or 2 = initially, while dark energy is either in potential regime 
Wd ~ — 1 or approximately in the scaling regime Wd 
case where dark energy starts in scaling regime, figure 



b 2 = — ot(f) 1 rather quickly. Although in the cases where <fi 2 = acf) 1 and <f> 2 



wj at initial stage. In the 
shows that <fi 2 converges to 
initially, 

4> 2 is subdominant and dark energy can be approximately in the scaling regime during 
the initial stage, the relation 2 ~ — a<p l is not satified. Hence, in these cases, SS is 
not constant during the initial stage as shown by thick long dashed and dashed lines in 
figure ([3]). However, after the initial transient, SS becomes constant. We also consider 



the case where Wd ~ — 1 and 



— a< 



initially. In this case SS is also constant in the 



early epoch. The above results for a^O also hold in the case of quintessence model. 



0.5 



O 



-0.5 




Figure 1. The evolution of Wd and the equation of state parameter w 2 of </> 2 for 
various initial conditions in the case of quintessence model. The tick lines represent 
Wd, while the thin lines represent W2- In this plot, Ai = 10, A2 = 1 and a = 0. 



From the above consideration, it can be concluded that in the early epoch if <fi 2 is 
subdominant compared with 1 , the entropy perturbations in dark energy SS tend to 
be constant on large scales, and consequently can survive until the present epoch for 
generic evolution of dark energy. The constancy of the entropy perturbations is not 
due to the particular form of Lagrangian ffTO]) . but is a consequence of |H| <C H and 
/jg <C H 2 , which can generally occur in two-field system when one field is subdominant 
compared with another. 



3.4- The effects of entropy perturbations on the CMB power spectrum 

Having consider the evolution of 5S in the early epoch,let us now study the evolution and 
effects of 5S about the present era. Since H and /i| depend on P s and P t xs, the evolution 
of SS and the influence of SS on Sd are governed by these quantities. However, since the 
contributions from P5 and P t xs are approximately the same, sometime we consider Pg 
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Figure 2. The evolution of wa and (ft 2 /(ft 1 for the case of k-essence model. The thick 
line represents the evolution of Wd, while the thin lines correspond to (ft 2 /(ft 1 . The tin 
solid, thin long dashed and thin dashed lines correspond to (ft 2 = —a^ 1 , (ft 2 — and 
(ft 2 = atft 1 at initial time respectively. In this plot, Ai = 50, A2 = 1 and a = 0.5. 
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Figure 3. The evolution of the long wavelength mode of 5S for models in figures 
([1]) and ([2]). The thin and thick lines correspond to the cases of a = and a = 0.5 
respectively. The thin long dashed line is the case where W2 ~ 1 while ss — 1 initially. 
The thick solid, thick long dashed and thick dashed lines correspond to (ft = —ot(ft\ 
(ft — and (ft 2 — aij) 1 at initial time respectively. 



only. As have been shown, \P t s/ ' Pd\ "C 1 in the early era. Using similar approach one 
expects that \P t s/Pd\ < 1 in the present epoch. This is because when 2 is dominant 
compared with 1 at present, we have |e 2 | > [e x | , | s 1 1 > |s 2 | and IP2I > |-P,i| so that 
\P t s\ ~ I-P.is 1 + P,2S 2 \ < |-P,2e 2 | ~ I Po- 1 ~ Pd- The ratio \P t s/pd\ can be of order of unity 
during the transition from 1 - dominant to 2 - dominant, because |s 2 | and IP2I are of 
order of | s 1 1 and |Pi| respectively. This rough estimate is confirmed by the numerical 
results in figure (jl]). Hence, 5S can have a contribution to the dynamics of 5d only during 
the transition stage, because the contribution from the entropy perturbations to the 
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pressure perturbations in (T4"T|) depends on the magnitude of Pg and Pxs- This implies 
that in the general two-field dark energy models, although the entropy perturbations in 
dark energy can be constant until the present epoch because one field is subdominant 
compared with the other in the early era, the entropy perturbations cannot significantly 
influence the dynamics of density perturbations in the universe unless the subdominant 
field becomes more important at late time. 

The evolution of 5S can influence the evolution of Sd via the relation 8PJ pd oc P t s8S 
given in flU]). Here, P 5 = [P S (1 + 4) ~ 2c 2 e XP iXS ] / p d . We follow CMBEASY code 
to normalize all perturbation variables such that TZ = — 1. If the value of P t s^S is 
sufficiently large and possitive, 5S will give a possitive contribution to SP, and also iid 
in f )34|) . Consequently, it follows from f l33j) that 5d can decrease more rapidly on large 
scales when P t s$S increases. This qualitative consideration is independent of the form 
of P(X, cj) 1 ) and in agreement with the numerical results in figure (jSJ). In contrast, if 
P t sSS is sufficiently large and negative, its contributions should increase 5j on large 
scales. 

For illustration, in figures @ and (jSJ), we plot the evolutions of P : s/p d , P,s$S and S d 
from the numerical integration for k-essence model whose Lagrangian is given by ( 12~TT) . 
Since the evolutions of these quantities for our choice of quintessence and k-essence 
models have approximately the same features, we consider the case of k-essence model 
only. 

We set Ai = 50 and A2 = 1 in the numerical integration. The initial conditions 
for the background field are chosen such that dark energy is in the scaling regime while 
2 = — a<p l initially. For the perturbed quantities, the adiabatic initial conditions given 
in ( 1541) are used for 5d and u d - The initial value of SS are chosen such that the CMB 
power spectrum at i = 2 for a = is maximally suppressed. This initial value can 
be in agreement with the initial value from inflation if \P t x\ < 1 during inflation. In 
our case, initial SS/TZ < corresponds to P } s5S/1Z < at the initial time, and leads 
to P^SS/TZ < during the transition stage. In general, the relation between the sign 
of P t s8S during the transition stage and the relative signs of 5S or P^s^S to 7Z at the 
initial time depends on the form of Lagrangian for dark energy. 

Since the effects of entropy perturbations in dark energy are most important on 
large scales we consider the perturbation mode whose wavenumber is k ~ 2 x 10~ 4 
Mpc -1 . This mode enters the Hubble radius about the present epoch. 

The effects of entropy perturbations on the density perturbations in the universe 
can be studied using the relation between the metric perturbations and 5d in f l25|) . 
Since we are interested in the perturbations on large scales and radiation contributes a 
negligible fraction of the total energy density during the transition stage, we write (125]) 
as 

$ = -H$ - y (n m 5 m + n d 5 d ) , (84) 

where the subscript m denotes matter. The relation $ ~ ^ has been used in the 
above equation because radiation is negligible. Hence, if the contributions from density 
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Figure 4. The evolution of P t s/pd- for various a. The solid, long dashed and dashed 
lines represent a = 0,0.5 and 0.8 respectively. 

15 r 
lO 

5 - 



O - 
-5 - 

_10 I ' ' ■ ' ' 1 

-3 -2.5 -2 -1.5 -1 -0.5 O 

log a 

Figure 5. The evolution of P t s5S (thick lines) and Sd (thin lines), the solid, long 
dashed and dashed lines correspond to a = 0, 0.5 and 0.8 respectively. For comparison 
, the evolution of Sd in the case where 5S = is plotted using dotted line. In this plot, 
k ~ 2 x 10~ 4 Mpc" 1 . 

perturbations is ignored for a while, it can be read from ( 184)) that the amplitude of metric 
perturbations decreases with time due to the expansion of the universe. The decreases 
of |$| and |^| which correspond to the ISW effect, partially cancel the ordinary Sachs- 
Wolfe (OSW) effect. It follows from ( )84|) that the inclusion of density perturbations can 
modify the evolution of $, such that if §/5d > 0, |$| will decrease faster, and |$| will 
decrease more slowly or increase if < [221 [231 [21]. When |$| decreases faster, 
the ISW contributions will be larger, leading to more cancellations between ISW and 
OSW effects. In this case, the CMB power spectrum is suppressed at low multipoles. In 
contrast, OSW effect is less canceled by ISW effect when |$| decreases more slowly, and 
OSW and ISW can add up when |$| increases. Consequently, the CMB power spectrum 
is enhanced at low multipoles in these cases. 



l 
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Using the same model of k-essence and initial conditions as above, we plot the 
evolutions of fld$d and $ in figure ([6]) and plot the CMB power spectrum in figure 
(J7j). In our numerical integration, $ is negative because all perturbation variables are 
normalized such that 1Z = — 1. Thus, in figure ([6]), $ increases when Qd$d < 0- This 
means that due to the contribution from 5d, \<p\ decreases more rapidly, i.e., the ISW 
effect is enhanced. The enhancement of ISW effect can lead to the suppression of the 
CMB power spectrum at low multipoles as shown in figure ([7]). These are in agreement 
with the above qualitative consideration 

1.0 , — , , , , , — , 




log a 



Figure 6. The evolution of fldSd is shown in the upper panel, while the evolution of 
<f> is shown in the lower panel. In the both panels, the solid, long dashed and dashed 
lines correspond to a — 0, 0.5 and 0.8 respectively. The evolution of both quantities 
in the case where SS = is denoted by dotted line. In this plot, k ~ 2 x 10~ 4 Mpc -1 . 

From the above analysis and figures (jHJ) - ©, one can conclude that if the curvature 
perturbation TZ and the contribution from entropy perturbations in dark energy, which 
is proportional to P : s^S, have opposite signs after the matter era, the ISW effect can 
be enhanced and consequently the CMB power spectrum can be suppressed at low 
multipoles. In contrast, if TZ and P } s5S have the same signs after the matter era, 
the CMB power spectrum at low multipoles can be enhanced. This is in agreement 
with literature (22J E3j. However, for two-field dark energy models, the contribution 
from entropy perturbations in dark energy depends on P t g, i.e., the evolution of dark 
energy, during the transition from <ft l - dominant to <f) 2 - dominant. Hence, in general, 
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Figure 7. The CMB power spectra for the k-essence model. The solid, long dashed 
and dashed lines correspond to a = 0, 0.5 and 0.8 respectively. The case where 5S = 
while a = is denoted by dotted line. 

P^SS/IZ < at the initial time does not necessarily lead to P^SS/IZ < after the 
matter era, and consequently the suppression of the CMB power spectrum at low 
multipoles. 

In order to study the influence of a on the effect of entropy perturbations, we plot 
the fractional difference in the CMB power spectra for k-essence and quintessence in 
figure (jHJ). The fractional difference in the CMB power spectra between the case of 
SS = and the case of SS ^ is expressed as the ratio 

AQ=^3, (85) 

where G\ is the power spectrum for the case of SS = and Ce is the power spectrum 
for the case of SS ^ 0. 

In this figure, the initial value for SS in the case of quintessence is also chosen 
such that the CMB power spectrum at £ = 2 for a = is maximally suppressed. 
This initial value of SS is much larger than a possible initial value from inflation given 
in (I66I) . because Px = 1 for quintessence. If the initial value in fl66|) is used, the 
entropy perturbations have small contribution to ISW effect and hence have no effect 
on CMB power spectrum. This figure shows that the suppression of the CMB power 
spectrum depends on a in the case of k-essence, but weakly depends on a in the case of 
quintessence. The influence of a on the suppression of CMB power spectrum strongly 
depends on the form of g(Y) in (jTUl) . 

4. Conclusion 

In this work, we study dynamics of entropy perturbations in two-field assisted dark 
energy model. Based on the scenario for this dark energy model in which one field 
is subdominant compared with the other in the early era, we show that the entropy 
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Figure 8. The fractional difference in the CMB power spectra for k-essence and 
quintessence models. 



perturbations in this two-field system tend to be constant on large scales in the early 
era, and hence can survive until the present epoch, for a generic evolution of dark energy 
during the radiation and matter eras. This behaviour of the entropy perturbations 
is preserved even when the fields are coupled via kinetic interaction. Although in 
our case the constancy of the entropy perturbations is obtained using the Lagrangian 
which has scaling behaviour, this result can also be obtained from more general form 
of Lagrangian. This is because the constancy of the entropy perturbations is the 
consequence of the subdomination of one field compared with the other in the early 
epoch. However, the entropy perturbations cannot significantly influence the dynamics 
of density perturbations in the universe unless the subdominant field in the early epoch 
becomes more important at late time. 

Since, for assisted dark energy, the subdominant field in the early epoch becomes 
dominant at late time, the entropy perturbations can influence the dynamics of metric 
perturbations so that the ISW effect is modified. The influences of the entropy 
perturbations on the ISW effect depend on both the initial condition for SS and 
the evolution of dark energy during the transition between 1 - domination and <p 2 - 
domination after the matter era. Assuming correlations between the entropy and 
curvature perturbations, the ISW effect will be enhanced and therefore the CMB 
power spectrum will be suppressed at low multipoles if the contribution from entropy 
perturbations and the curvature perturbations have opposite signs, i.e., P^SS/TZ < 0, 
after the matter era. For our choice of models in the numerical integration, the negative 
P t sSS/lZ at initial time leads to negative P^SS/1Z after the matter era. However, in 
general, initial P^SS/TZ < does not necessarily imply P^SS/TZ < after the matter 
era. 

Assuming dark energy to exist since inflationary epoch, the initial value of \SS\ 
estimated during inflationary epoch cannot be large enough to influence the ISW effects 
in the case of quintessence. This is because P >x = 1 for usual quintessence models. The 
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initial value of \SS\ can be large enough to influence the ISW effect if \P t x\ is significantly 
smaller than unity during inflationary era, which is possible in many k-essence models. 

If dark energy fields are coupled via kinetic terms, the modifications of CMB power 
spectrum due to entropy perturbations also depend on the coupling strength a. For our 
choice of models in the numerical integration, this dependence is significant in the case 
of k-essence, but not in the case of quintessence. The influence of a on the CMB power 
spectrum strongly depends on the form of g(Y) in (TTOT) . 

Acknowledgments 

The author would like to thank A. Ungkitchanukit and anonymous referee for comments 
on the manuscript. This work is supported by Thailand Research Fund (TRF) through 
grant TRG5280036. 

[I] A. G. Riess et al. [Supernova Search Team Collaboration], 2004 Astrophys. J. 607 665 

[arXiv : astro-ph/0402512] ; 

R. A. Knop et al., [Supernova Cosmology Project Collaboration], 2003 Astrophys. J. 598 102 

[arXiv : astro-ph/0309368] . 
[2] C. L. Bennett et al., 2003 Astrophys. J. Suppl. 148 1 [arXiv: astro-ph/0302207]. 
[3] M. Tegmarket al. [SDSS Collaboration], 2004 Phys. Rev. D69 103501 [arXiv : astro-ph/0310723]. 
[4] J. Frieman, M. Turner and D. Huterer, 2008 Ann. Rev. Astron. Astrophys. 46 385 

[arXiv: 0803. 0982]. 

[5] C. W. Chen, Je-An Gu and P. Chen, 2009 Mod. Phys. Lett. A24 1649 [arXiv: 0903. 2423]. 

[6] C. Wetterich, 1988 Nucl. Phys. B302 668. 

[7] B. Ratra and J. Peebles, 1988 Phys. Rev. D37 321. 

[8] R. R. Caldwell, R. Dave and P. J. Stcinhardt, 1998 Phys. Rev. Lett. 80 1582 

[arXiv : astro-ph/9708069] . 
[9] T. Chiba, T. Okabe and M. Yamaguchi, 2000 Phys. Rev. D62 023511 [arXiv : astro-ph/9912463]. 
[10] C. Armendariz-Picon, V. Mukhanov, and P. J. Steinhardt, 2000 Phys. Rev. Lett. 85 4438 

[arXiv : astro-ph/0004134] . 

[II] P. Kanti and K. A. Olive, 1999 Phys. Rev. D60 043502 [arXiv:hep-ph/9903524]. 

[12] P. G. Ferreira and M. Joyce, 1997 Phys. Rev. Lett. 79 4740 [arXiv: astro-ph/9707286]. 

[13] E. J. Copeland, A. R. Liddle and D. Wands, 1998 Phys. Rev. D57 4686 [arXiv:gr-qc/9711068]. 

[14] A. R. Liddle, A. Mazumdar and F. E. Schunck, 1998 Phys. Rev. D58 061301 

[arXiv : astro-ph/9804177] . 
[15] S. A. Kim, A. R. Liddle and S. Tsujikawa, 2005 Phys. Rev. D72 043506 

[arXiv : astro-ph/0506076] . 
[16] J. Ohashi and S. Tsujikawa, 2009 Phys. Rev. D80 103513 [arXiv: 0909. 3924]. 
[17] D. A. Easson, R. Gregory, D. F. Mota, G. Tasinato and I. Zavala, 2008 JCAP 0802 010 

[arXiv: 0709. 2666]. 

[18] L. P. Chimento, M. Forte, R. Lazkoz and M. G. Richarte, 2009 Phys. Rev. D79 043502 
[arXiv: 081 1.3643]. 

[19] C. van de Bruck and J. M. Weller, 2009 Phys. Rev. D80 123014 [arXiv: 0910. 1934]. 
[20] E. N. Saridakis and J. M. Weller, 2010 Phys. Rev. D81 123523 [arXiv : 0912 . 5304]. 
[21] E. N. Saridakis and S. V. Sushkov, 2010 Phys. Rev. D81 083510 [arXiv: 1002.3478]. 
[22] C. Gordon and W. Hu, 2004 Phys. Rev. D70 083003 [astro-ph/0406496]. 
[23] T. Moroi and T. Takahashi, 2004 Phys. Rev. Lett. 92 091301 [arXiv: astro-ph/0308208]. 
[24] K. Karwan, 2007 JCAP 0707 009 [astro-ph/0701009]. 

[25] K. A. Malik and D. Wands, 2005 JCAP 0502 007 [arXiv: astro-ph/0411703]. 



Dynamics of entropy perturbations in assisted dark energy with mixed kinetic terms 25 

[26] C. Gordon, D. Wands, B. A. Bassett and R. Maartens, 2000 Phys. Rev. D63 023506 

[arXiv : astro-ph/0009131] . 
[27] D. Langlois and S. Renaux-Petel, 2008 JCAP 0804 017 [arXiv: 0801 . 1085]. 
[28] D. Langlois, S. Renaux-Petel, D. A. Steer and T. Tanaka, 2008 Phys. Rev. Lett. 101 061301 

[arXiv: 0804. 3139]. 

[29] D. Langlois, S. Renaux-Petel, D. A. Steer and T. Tanaka, 2008 Phys. Rev. D78 063523 
[arXiv: 0806. 0336]. 

[30] R. Bean, S. H. Hansen and A. Melchiorri, 2001 Phys. Rev. D64 103508 [arXiv : astro-ph/0104162]. 
[31] F. Piazza and S. Tsujikawa, 2004 JCAP 0407 004 [arXiv : hep-th/0405054] . 
[32] J. M. Bardeen, 1980 Phys. Rev. D22 1882. 

[33] A. J. Christophcrson and K. A. Malik, 2009 Phys. Lett. B675 159 [arXiv: 0809. 3518]. 
[34] A. De Felice, S. Mukohyama and S. Tsujikawa, 2010 [arXiv: 1006.0281]. 

[35] M. Doran and C. M. Miiller, G. Schafer and C. Wetterich, 2003 Phys. Rev. D68 063505 

[astro-ph/0304212] . 
[36] K. A. Malik, 2001 [arXiv: astro-ph/0101563]. 

[37] A. Riotto, 2002 Trieste 2002, Astroparticle physics and cosmology* 317 [hep-ph/0210162]. 
[38] M. Doran, 2005 JCAP 0510 011 [arXiv: astro-ph/0302138]. 



